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ABSTRACT 


This  research  focuses  on  the  detection  of  a  low  power  communication  signal  in  the 
presence  of  a  strong  co-channel  television  broadcast  interference  signal.  The  presence  of 
strong  co-channel  interference  makes  the  recovery  of  the  desired  weak  power  signal 
impossible  using  conventional  filtering  techniques  that  are  based  on  time  and  frequency 
characteristics  of  the  signals.  The  second-generation  digital  video  broadcasting  terrestrial 
(DVB-T2)  standard  is  employed  as  co-channel  interference  in  an  additive  white  Gaussian 
noise  channel.  The  weak  signal  is  assumed  to  have  a  considerably  smaller  bandwidth  than 
the  TV  interference  and  negligible  phase-shift  due  to  multipath. 

By  using  two  antennas  at  the  receiver,  channel  diversity  can  be  exploited,  and  the 
weak  signal  can  be  recovered  using  a  Kalman  filter  (KF),  assuming  the  channels  seen  by 
the  two  antennas  are  independent  and  time-invariant.  Moreover,  the  transmitted  co¬ 
channel  interference  is  modeled  as  the  state  of  a  dynamic  system  whose  input  is  the 
signal  received  at  one  antenna  and  the  output  is  the  signal  received  at  the  second  antenna. 
Within  this  framework,  the  state  can  be  estimated  by  a  KF.  Channel  estimation  is 
performed  using  DVB-T2  pilots.  Performance  of  the  system  is  evaluated  at  different 
signal-to-noise  ratio  (SNR)  and  signal-to-interference  ratio  (SIR).  Results  show  that  the 
weak  signal  can  be  reconstructed  with  bit  error  ratio  (BER)  of  10-^  or  less  under  most 
SNR  and  SIR  conditions  considered  in  the  study. 
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EXECUTIVE  SUMMARY 


Most  conventional  signal  processing  techniques  used  in  receivers  focus  on  signal  time 
and  frequency  characteristics  to  recover  transmitted  information.  In  such  cases,  signal 
information  can  only  be  recovered  if  its  power  is  greater  than  the  combined  power  of  the 
noise  and  interferences;  however,  military  and  commercial  application  scenarios  exist  in 
which  the  desired  signal  is  transmitted  in  the  presence  of  much  stronger  co-channel 
interference  signal. 

The  objective  of  this  thesis  is  to  recover  a  signal  embedded  in  strong  co-channel 
broadcasting  interference  signal.  The  overhead  information  of  the  broadcasting  signal 
(typically  used  for  detection  and  synchronization  tasks)  is  used  in  a  multi-antenna  system 
to  estimate  the  multipath  and  time  delay  characteristics,  allowing  for  a  relatively  simple 
interference  cancellation  algorithm  based  on  Kalman  filtering  (KF)  techniques.  The 
technique  presented  is  novel  and  may  open  the  possibility  of  using  the  frequency 
spectrum  occupied  by  digital  video  and  radio  broadcasting  for  other  purposes,  such  as 
low  power  data  transmission. 

The  approach  considered  in  this  work  is  based  on  a  two-antenna  model  at  the 
receiver  and  relies  on  channel  diversity  at  the  two  receiving  antennas.  The  main 
assumptions  of  this  approach  are: 

1 .  The  weak  signal  is  minimally  affected  by  multipath. 

2.  The  knowledge  of  the  interferer’s  overhead  is  sufficient  to  estimate  the 
transmission  channel  parameters. 

These  conditions  can  be  satisfied  by  placing  the  two  antennas  at  least  half  a  wavelength 
apart  and  keeping  the  weak  signal  bandwidth  smaller  than  the  frequency  coherence  of  the 
channel. 

In  this  research,  the  strong  signal  interferer  is  assumed  to  be  a  second-generation 
digital  video  broadcasting-terrestrial  (DVB-T2)  signal,  which  is  based  on  orthogonal 
frequency-division  multiplexing  (OFDM),  and  the  weak  signal  is  a  quadrature  phase-shift 
keying  (QPSK)  single-carrier  modulated  signal.  The  DVB-T2  signal  passes  through  two 
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multipath  channels  with  additive  white  Gaussian  noise  (AWGN).  Estimation  of  the  two 
channels  is  conducted  using  the  pilots  present  in  the  DVB-T2  signal. 

The  proposed  approach  is  based  on  the  KF  and  uses  the  signals  received  at  the 
two  antennas  and  estimated  channel  coefficients  to  estimate  the  DVB-T2  signal 
interferer.  By  this  particular  structure,  adapted  from  the  controls  literature,  the  two 
signals  received  at  the  antennas,  respectively,  act  as  predictor  and  corrector  of  the  state  of 
a  dynamic  system,  where  the  state  represents  the  transmitted  TV  broadcasted  signal. 
What  the  KF  literature  calls  the  innovation,  contains  the  unpredictable  part  of  the 
dynamic  model,  from  which  the  weak  signal  can  be  recovered  and  demodulated  using  a 
standard  QPSK  demodulator.  This  concept  is  illustrated  in  Figure  1 . 


Figure  1.  Block  diagram  of  the  KF  approach  used  to  recover  the  weak  signal 


The  estimated  channel  impulse  responses  (CIRs)  are  used  both  for  equalization  in 
the  respective  DVB-T2  receivers  and  in  the  KF  to  estimate  the  state  of  the  DVB-T2 
signal. 
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Results  show  that  the  weak  signal  can  be  extracted  from  the  interference  and 
noise  when  the  underlying  assumptions  of  the  scheme  presented  in  the  research  are  met. 
Both  interference  and  weak  signals  can  be  demodulated  with  bit  error  ratio  (BER)  values 
of  10-  or  less  under  most  signal-to-noise  ratio  (SNR)  and  signal-to-interference  ratio 
(SIR)  conditions  investigated.  Follow-on  research  should  include  a  robust  approach  for 
estimating  the  channel  characteristics  to  allow  for  better  BER  performances. 
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I.  INTRODUCTION 


A.  BACKGROUND 

In  most  communications  systems,  conventional  signal  proeessing  techniques  used 
in  reeeivers  focus  on  time  and  frequency  characteristics  of  the  signals  to  reeover  the 
transmitted  information.  In  such  cases,  the  signal  ean  only  be  reeovered  if  its  power  is 
sufficiently  greater  than  that  of  the  noise  and  interferenee;  however,  there  may  be  eases 
when  the  desired  signal  power  is  eonsiderably  lower  (whether  intentionally  or 
unintentionally)  and  is  transmitted  in  the  presence  of  strong  co-channel  interference.  Sueh 
situations  can  be  encountered  in  both  military  and  commercial  systems.  For  example, 
military  communication  systems  may  experienee  strong  eo-ehannel  jamming,  whieh 
denies  the  reception  of  a  desired  signal.  Similarly,  in  the  commereial  world,  if  a 
communieation  system  ean  operate  reliably  in  the  presence  of  a  strong,  eo-ehannel 
interferenee,  then  there  is  no  need  to  buy  spectrum,  an  expensive  commodity. 

In  this  thesis  we  address  the  problem  a  reeovering  a  communication  signal  in  the 
presence  of  a  strong  co-ehannel  interferer.  In  partieular,  we  show  that  using  a  priori 
knowledge  of  the  interferer  (e.g.,  its  overhead  and  a  multi-antenna  reeeiver)  the  interferer 
can  be  estimated  and  highly  attenuated  by  eancellation. 

B.  OBJECTIVE 

The  objeetive  of  this  thesis  is  to  recover  a  signal  from  a  eo-channel  interference 
signal.  In  this  ease,  classical  filtering  techniques  eannot  be  applied,  as  both  signals 
occupy  the  same  frequency  bands.  The  signal-of-interest  is  extracted  by  taking  advantage 
of  the  known  overhead  portion  (eommonly  used  for  detection  and  synchronization) 
contained  in  the  interferer,  selected  to  be  a  broadcast  signal.  This  overhead  information 
can  be  used  in  a  multi-antenna  system  to  estimate  the  transmission  ehannel  multipath  and 
time  delay  eharacteristics,  whieh  leads  to  a  relatively  simple  interference  cancellation 
seheme  based  on  Kalman  filtering  teehniques. 
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The  technique  presented  in  this  work  is  novel  and  may  open  the  possibility  of 
using  the  frequency  spectrum  occupied  by  digital  video  and  radio  broadcasting  for  other 
purposes,  such  as  low  power  data  transmission. 

C.  RELATED  WORK 

The  co-channel  interference  problem  has  been  investigated  recently  for  various 
kinds  of  settings.  One  approach  has  been  based  on  Independent  Component  Analysis 
(ICA),  where  the  signals’  statistical  properties  are  used  for  separation  purposes  [1]. 
Algorithms  based  on  ICA  were  presented  in  [2]  which  show  the  effectiveness  of  the 
technique  in  cases  with  no  multipath. 

The  use  of  Kalman  Filtering  in  a  multi-antenna  system  has  recently  been 
introduced  in  [3]  and  [4].  The  use  of  the  Kalman  Filter  (KF)  for  this  particular  application 
has  not  been  widely  published  and  seems  to  yield  an  interesting  framework  for 
recovering  data  from  a  multi-antenna  system.  Formulating  the  multi-antenna  receiver  as  a 
state  space  dynamic  system,  widely  used  in  the  controls  systems  literature,  allows  for  the 
application  of  well-known  state  estimation  techniques,  such  as  the  state  observer  and  the 
KF. 

D.  APPROACH 

The  approach  employed  in  this  research  is  based  on  a  two-antenna  model  at  the 
receiver  and  relies  on  the  presence  of  channel  diversity  at  the  two  receiving  antennas. 
This  means  that  the  signal  reaching  each  antenna  sees  a  different  channel.  The  main 
assumptions  of  this  approach  are  the  following: 

1 .  The  weak  signal  is  minimally  affected  by  the  multipath. 

2.  Knowledge  of  the  interferer’s  overhead  allows  for  the  estimation  of  the 
transmission  channel  parameters. 

These  conditions  can  be  satisfied  by  placing  the  two  antennas  at  least  half  a 
wavelength  apart  and  keeping  the  weak  signal  bandwidth  smaller  than  the  frequency 
coherence  of  the  channel. 
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The  strong  signal  interferer  is  assumed  to  be  a  second  generation  digital  video 
broadcasting-terrestrial  (DVB-T2)  signal,  which  is  the  European  television  (TV) 
standard.  DVB-T2  is  based  on  orthogonal  frequency-division  multiplexing  (OFDM).  The 
weak  signal  is  modeled  as  a  quadrature  phase-shift  keying  (QPSK),  single-carrier 
modulated  signal.  The  DVB-T2  signal  passes  through  two  additive  white  Gaussian  noise 
(AWGN)  channels.  Estimation  of  the  two  channels  characteristics  is  conducted  via  the 
pilots  present  in  the  DVB-T2  signal. 

The  proposed  approach  is  based  on  the  KF  algorithm  and  uses  the  signals 
received  at  the  two  antennas  and  estimated  channels  characteristics  to  estimate  the  DVB- 
T2  signal  interferer.  By  this  particular  structure,  adapted  from  the  controls  literature,  the 
two  signals  received  at  the  antennas  act  as  predictor  and  corrector  of  the  state  of  a 
dynamic  system,  where  the  state  represents  the  transmitted  TV  broadcasted  signal.  What 
is  called  the  “innovation”  in  the  KF  literature  contains  what  cannot  be  predicted  by  the 
dynamic  model,  from  which  the  weak  signal  can  be  recovered  and  demodulated  using  a 
standard  QPSK  demodulator.  This  concept  is  illustrated  in  Figure  1.  The  performance  of 
the  system  is  evaluated  at  different  signal-to-noise  ratio  (SNR)  and  signal-to-interference 
ratio  (SIR)  values. 

E.  THESIS  ORGANIZATION 

This  thesis  is  organized  into  five  chapters  including  the  introduction.  The  basic 
idea  behind  the  DVB-T2  structure  is  explained  in  Chapter  II  to  highlight  the  estimation 
process  of  the  channel  impulse  response  (CIR).  The  KF  algorithm  and  its  application  to 
the  proposed  scheme  are  presented  in  Chapter  III.  The  discussion  in  Chapter  IV  focuses 
on  the  simulation  model  and  results  obtained  at  different  SNRs  and  SIRs.  Finally, 
conclusions  and  recommendations  are  provided  in  Chapter  V. 
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II.  DVB-T2  STRUCTURE 


We  have  seen  the  evolution  TV  has  undergone  over  the  last  couple  of  decades. 
The  video  which  we  see  on  standard  definition  television  (SDTV)  follows  standards  that 
were  developed  in  the  1940s  and  1950s.  The  two  most  widely  used  analog  standards  are 
the  national  television  systems  committee  (NTSC)  and  phase  alternating  line  (PAL). 
NTSC  is  the  dominant  standard  on  the  American  continent  and  Japan,  while  PAL  is  the 
one  for  Europe,  Asia,  and  Australia.  These  two  standards  have  dominated  the  TV  market 
for  almost  60  years.  With  the  advent  of  the  digital  age,  NTSC  and  PAL  have  lost 
considerable  ground  to  emerging  technologies. 

Digital  TV  (DTV)  broadcasting  concepts  are  presented  in  this  chapter  followed  by 
a  brief  description  of  OFDM  and  its  implementation  in  the  DVB-T2  standard.  Finally,  we 
describe  the  approach  followed  to  extract  channel  characteristics  using  DVB-T2  pilot 
signals. 

A.  DIGITAL  TELEVISION 

Technology  advances  have  taken  huge  leaps  in  the  recent  past,  and  tremendous 
progress  has  been  seen  in  every  field  of  life.  Such  advances  have  resulted  in  rising 
demand  for  fast,  efficient,  and  high  throughput  systems,  and  the  broadcasting  industry  is 
no  exception.  In  order  to  keep  up  with  the  technology  and  consumer  demand,  the  shift  to 
DTV  is  a  natural  evolution.  The  first  DTV  service  was  delivered  via  satellite,  and  now 
this  service  is  provided  via  terrestrial  and  cable  networks.  DTV  has  many  advantages 
over  analog  TV.  A  Few  of  these  advantages  include: 

•  less  interference 

•  better  video  and  audio  quality 

•  multiplexing  which  allows  various  types  of  signals  to  be  transmitted 
simultaneously 

•  better  spectral  efficiency 
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Due  to  the  inherent  versatility  of  digital  transmission,  DTV  has  brought  many 
exciting  variations  to  TV  broadcasting,  and  various  digital  high-defmition  TV  (HDTV) 
standards  have  been  adopted  all  over  the  world.  The  most  popular  ones  are: 

•  advanced  Television  Systems  Committee  (ATSC) 

•  integrated  services  digital  broadcasting  -  terrestrial  (ISDB-T) 

•  digital  video  broadcasting  (DVB-T) 

•  digital  terrestrial  multimedia  broadcast  (DTMB) 

The  regions  where  these  standards  have  been  implemented  are  shown  in  Figure  2. 
All  these  standards  follow  the  same  general  structure;  however,  specific  implementation 
schemes  differ  from  standard  to  standard. 


Q/Qfl  ■  CXGra  ■  ATSC  ■  SOe-T  B  DTWR 

Figure  2.  Map  of  digital  terrestrial  television  broadcast  standards, 

from  [5]. 


The  research  conducted  here  considers  the  case  of  a  strong  TV  broadcast  signal 
interfering  with  a  signal  transmitting  at  a  lower  power  embedded  in  the  same  frequency 
band.  Although  the  assumption  throughout  is  that  the  TV  broadcast  signal  has  a  DVB-T2 
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structure,  the  techniques  developed  are  general  in  nature  and  can  be  applied  to  any 
standard. 

The  DVB-T2  standard  structure  is  based  on  OFDM  and  employs  a  number  of 
digital  modulation  constellations  such  as  QPSK,  16-quadrature  amplitude  modulation 
(QAM),  64-QAM,  or  256-QAM,  and  various  coding  schemes  to  adapt  to  different 
channel  and  multipath  situations. 

B.  ORTHOGONAL  FREQUENCY-DIVISION  MULTIPLEXING 

In  wireless  systems,  the  signal  at  the  receiver  comes  from  different  paths  due  to 
reflection  and  diffraction  from  large  objects  such  as  buildings,  hills,  etc.  Each  path  is 
characterized  by  different  levels  of  attenuation  and  delay;  therefore,  multiple  copies  of 
the  same  signal  are  received,  delayed,  and  attenuated,  causing  inter-symbol 
interference  (ISI). 

A  wireless  channel  is  characterized  by  its  time  spread  and  frequency  spread.  The 
effect  of  ISI  is  negligible  when  the  symbol  duration  of  the  transmitted  signal  is 
considerably  higher  than  the  time  spread  of  the  channel.  This  condition  corresponds  to  a 
flat  fading  channel  scenario;  however,  equalization  is  essential  to  compensate  for  channel 
distortions  when  the  transmitted  signal  symbol  duration  is  smaller  than  the  channel  time 
spread  value.  This  second  scenario  corresponds  to  a  frequency-selective  channel  case. 
Equalization  requires  training  data  (thus,  loss  of  bandwidth)  that  can  be  computationally 
expensive  if  done  blindly  and  is  difficult  when  the  channel  characteristics  are  time 
varying. 

In  single-carrier  modulation  (SCM),  each  bit  is  transmitted  on  a  single  carrier 
frequency.  An  alternative  approach  is  to  divide  the  available  channel  bandwidth  into  a 
number  of  sub-channels  each  with  equal  bandwidth.  The  bandwidth  of  each  channel  is 
selected  to  be  as  narrow  as  possible  so  that  each  sub-channel  sees  a  flat  fading  channel 
and  ISI  is  minimized.  Different  information  symbols  can  be  transmitted  on  each  sub¬ 
channel  simultaneously.  This  scenario  is  a  particular  case  of  frequency-division 
multiplexing  (FDM),  known  as  multicarrier  modulation  (MCM)  in  which  each  sub- 


7 


channel  is  associated  with  a  single  carrier.  SCM  and  MCM  concepts  are  illustrated  in 
Figure  3. 
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Figure  3.  SCM  and  MCM,  from  [6]. 


In  SCM,  data  is  transmitted  on  a  single  frequency  at  a  higher  data  rate  R  =  l/T. 
In  MCM,  the  data  is  transmitted  over  N  subcarriers  at  a  lower  data  rate  R/N  =  l/ NT, 
as  shown  in  Figure  4.  In  such  a  case,  2l  N-  symbol  data  block  can  be  chosen  so  that  the 
total  time  duration  NT  of  the  symbol  frame  is  considerably  longer  than  the  channel  time 
spread. 


Figure  4.  Block  diagram  of  an  MCM  modulation  scheme. 


8 


OFDM  is  a  special  case  of  FDM  in  which  all  subcarriers  are  orthogonal  to  each 
other  and  is  widely  used  in  wireless  systems  including  the  popular  IEEE  802.11a,  g,  n 
standards,  DVB-T2  standard  and  long-term  evolution  (LTE)  cellular  networks.  A  further 
advantage  of  MCM  is  that  the  transmission  power  can  be  efficiently  distributed  across 
subcarriers  based  on  the  transmission  channel  frequency  response  to  maximize  the  data 
rate. 

As  shown  in  Figure  4,  data  bits  transmitted  at  a  rate  R  =  \IT  bits  per  second  are 
converted  into  N  streams,  each  one  having  a  symbol  rate  of  F^=\l NT  =  \IT^ 
symbols  per  second  (sps)  and  independently  modulated  by  one  of  N  orthogonal 
subcarriers.  By  separating  the  subcarriers  by  MT^ ,  we  guarantee  orthogonality  so  that  the 

transmitted  symbols  in  each  time  data  block  do  not  interfere  with  each  other.  Note  that  it 
is  not  feasible  to  implement  the  scheme  given  in  Figure  4  due  to  the  large  number  of 
modulators  used.  Instead  baseband  implementation  is  considered,  which  is  discussed 
next  [7]. 

1.  Baseband  Implementation 

It  is  well  known  that  an  OFDM  signal  is  transmitted  as  a  sequence  of  OFDM 
frames,  each  one  of  duration  and  where  each  frame  .^(t)  is  the  superposition  of 
modulated  subcarriers,  as  given  by 

=  (2.1) 

^  k=Q 

where  (t)  is  given  by 

=  0<t<T,  (2.2) 

with  k  =  0,1,2,....A-1,  and  represents  the  complex  envelope  of  the  N’'  subcarrier, 

is  the  nominal  carrier  frequency,  and  iSf  -MT^  is  the  carrier  symbol  rate  chosen  to 

maintain  orthogonality  among  subcarriers.  The  complex  envelope  of  the  k*''  subcarrier  is 
given  by 

z^{t)^Z{k)^I^+jQ^,  0<t<r^, 
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(2.3) 


where  I,^  and  represent  /,  Q  values  of  the  k‘^  subcarrier  symbol  and  are  the  in-phase 
and  quadrature  components  of  (t),  respectively. 

By  substituting  Equation  (2.2)  into  Equation  (2.1),  we  get 

N  k=o  N  k=Q 

0<t<T^.  (2.4) 

Thus,  the  complex  envelope  of  ODFM  signal  is  given  by 

=  0<(<r,.  (2.5) 

^  k=0 


=  Re 


N 


JlTrkAft 


Sampling  rate  and  sampling  time  are  given  by  =  NF^  and  =TJN, 

which  yields  the  sampled  signal  expression: 

^  nTs^ 


V  N  j 


1  A-l 

=  -Yz(k)e' 

k=0 


JlTikAfnTJ  N 


n  =  0,l,....,7V-l. 


(2.6) 


Choosing  A/T^  =  1  renders  the  subcarriers  orthogonal  to  each  other  and  yields  the 
well-known  expression  of  the  inverse  discrete  Fourier  transform  (IDFT) 

jlTvknlN 


k=Q 

=  IDFT{Z(yS:)}. 


n  =  0,l,....,iV-l 


(2.7) 


It  is  evident  from  Equation  (2.7)  that  the  sequence  s^{n)  can  be  obtained  by 
taking  the  IDFT  of  Z{k^  where A:  =  0,I,....,A^-I.  Conversely,  Z{k^  can  be  obtained  by 
taking  the  discrete  Fourier  transform  (DFT)  of  the  sequence  Sj^{n),  as  given  by 

Z(A:)  =  ^^^  A:  =  0,1,. ...,7V -I 

k=0 

=  DFT{5^(n)}.  (2.8) 
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Both  DFT  and  IDFT  operations  can  be  computed  via  very  efficient  fast  Fourier 
transform  (FFT)  and  inverse  fast  Fourier  transform  (IFFT)  algorithms,  respectively. 
Therefore,  OFDM  modulation  and  demodulation  can  be  implemented  using  the  IFFT  and 
FFT  algorithms,  respectively,  at  baseband.  In  practice,  the  number  of  subcarriers  N  is 
chosen  so  that  the  FFT  operation  can  be  implemented  efficiently.  Although  in  most 
implementations  a  power  of  four  is  preferred  (all  IEEE  standards),  it  is  not  necessarily  the 
case  in  all  recent  implementations,  such  as  in  DTMB  adopted  by  the  People’s  Republic  of 
China  [8]. 

A  typical  implementation  of  an  OFDM  modulator  via  baseband  approach  is 
shown  in  Figure  5.  Incoming  data  bits  {0,1}  are  encoded  into  M-QAM  complex  symbols, 
which  are  multiplexed  into  N  parallel  streams.  For  the  bit  stream,  a  group  of  log^M 

bits  is  mapped  into  a  subcarrier  symbol  Z{k)  =  /^  +  jQj^,  0  <  t  .  Next,  the  sequence  of 
A  symbols  is  fed  into  the  IFFT  block  to  obtain  the  sequence  which  represents  the 
samples  of  the  OFDM  complex  envelope. 


-sin  Infj 


Figure  5.  Block  diagram  of  OFDM  modulator,  from  [5]. 


The  digital-to-analog  (D/A)  converter  output  is  a  complex  I-Q  waveform  given  by 
=  N),  (2.9) 

/t=0 


where  p(t)  is  the  interpolation  pulse  shape,  coming  from  the  Zero  Order  Hold  and  the  low 
pass  reconstruction  filter  [5]. 
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2. 


Guard  Interval 


A  guard  interval  (GI)  is  added  to  an  OFDM  symbol  in  order  to  avoid  ISI.  The  GI 
duration  is  chosen  so  that  where  is  the  maximum  time  spread  of  the 

channel.  At  the  receiver,  GI  is  removed  before  the  FFT  operation.  Adding  GI  decreases 
the  throughput  of  the  system  because  no  new  information  is  present  in  the  GI.  The  effect 
of  r  and  GI  on  ISI  is  shown  in  Figure  6.  It  can  be  seen  that  ISI  cannot  be  eliminated 
when  <T . 
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Figure  6.  Effect  of  T  and  GI  on  ISI. 

Various  OFDM-based  techniques  have  been  established  based  on  the  structure  of 
the  data  used  in  GI.  Among  them,  Cyclic  prefix  (CP-OFDM),  zero  prefix  (ZP-OFDM), 
known  symbol  padding  (KSP-OFDM)  and  time-domain  synchronous  (TDS-OFDM)  are 
the  most  widely  used.  Next,  we  discuss  CP-OFDM  as  it  is  used  in  DVBT-2. 

3.  CP-OFDM 

In  this  technique,  the  GI  data  is  created  by  repeating  the  last  L  entries  of  each 

OFDM  frame,  so  that  LT^  exceeds  the  time  spread  of  the  channel,  as  shown  in  Figure  7. 
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This  constraint  forces  periodic  boundary  conditions,  and  it  can  be  easily  shown  that  each 
received  OFDM  frame  is  the  result  of  a  circular  convolution  between  the  transmitted 
frame  and  CIR,  thus  making  the  use  of  the  FFT  appropriate  [5]. 


Next,  demodulation  of  each  OFDM  frame  is  obtained  by  discarding  the  first  L 
samples,  eliminating  inter  block  interference,  and  then  taking  the  A^-point  DFT  of  the 
remaining  samples.  This  operation  yields 

Y{k)  =  H{k)X{k),  k  =  Q,l,-,N-l,  (2.10) 

where  H{k)  is  the  N  -point  DFT  of  the  CIR.  Therefore,  the  CIR  can  be  determined  by 
transmitting  training  data  at  the  beginning  of  the  transmission  and  tracked  by  reserving 
some  subcarriers  as  pilots  so  that  the  transmitted  data  ^  can  be  estimated. 

The  approach  considered  in  this  work  to  separate  the  desired  signal  from  the  TV 
broadcasting  interferer  makes  full  use  of  the  information  contained  in  the  DVB-T2 
signal’s  overhead.  Such  overhead  information  is  readily  available  and  easy  to  demodulate 
since  broadcasting  systems  are  designed  for  mass  production  and  require  quick 
synchronization  at  the  receivers. 

In  this  section,  we  present  the  characteristics  of  the  DVB-T2  system  with 

particular  attention  to  the  overhead,  which  contains  the  necessary  information  for 
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synchronization,  channel  estimation,  and  demodulation.  As  is  seen  in  following  chapters, 
the  proposed  system  is  very  efficient  in  the  sense  that  it  makes  minimal  use  of  the 
overhead  without  the  need  to  demodulate  the  complete  broadcast  signal. 


4.  DVB-T2 

DVB-T2  uses  coded  OFDM  (COFDM)  with  a  large  number  of  subcarriers  and 
employs  different  modulation  schemes  along  with  inner  low  density  parity  check  (LDPC) 
coding  combined  with  outer  Bose-Chaudhuri-Hocquengham  (BCH)  coding.  The  number 
of  carriers,  GI  sizes,  and  pilot  signals  can  be  adjusted.  Some  of  the  adjustable  parameters 
are  given  in  Table  1.  Full  details  of  the  DVB-T2  standard  can  be  found  in  [9]. 


Table  1 .  Parameters  used  in  DVB-T2  standard. 


FEC 

LDPC  +  BCH 

1/2,  3/5,  2/3,  3/4,  4/5,  5/6 

Modes 

QPSK,  16-QAM,  64-QAM,  256-QAM 

GI 

1/4,  19/128,  1/8,  19/256,  1/16,  1/32,  1/128 

FFT  size 

Ik,  2k,  4k,  8k,  16k,  32k 

Scatter  pilots 

1%,  2%,  4%,  8%  of  total 

Continual  pilots 

0.4%-2.4%  (0.4%  -  0.8%  in  8k-32k) 

Bandwidth 

1.7,  5,  6,  7,  8,  10  MHz 

5.  DVB-T2  Structure 

DVB-T2  frame  structure  is  shown  in  Figure  8.  At  the  top  level,  there  are 
superframes  that  consist  of  T2  frames,  which  are  further  divided  into  OFDM  symbols. 
The  number  of  T2  frames  in  a  super  frame  is  a  configurable  parameter. 

Preamble  PI  symbols  are  used  for  fast  detection  of  the  T2  signal  and 
identification  of  its  parameters  from  the  preamble  to  enable  proper  synchronization  and 
demodulation  at  the  receiver.  The  P2  symbol(s)  follow  immediately  after  the  PI  symbol. 
The  main  purpose  of  the  P2  symbol(s)  is  to  carry  LI  signaling  data  which  define  the 
various  parameters  used  for  the  transmission  [9]. 
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Each  T2  frame  comprises  of  one  PI  preamble  symbol,  followed  by  one  or  more 
P2  preamble  symbols,  followed  by  a  configurable  number  of  data  symbols. 
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Figure  8.  DVB-T2  frame  structure,  showing  the  division  into 
superframes,  T2  frames,  and  OFDM  symbols,  after  [7]. 


6.  DVB-T2  Pilots 

Various  cells  within  the  OFDM  frame  are  modulated  with  reference  information 
whose  transmitted  value  is  known  to  the  receiver.  These  are  called  pilots  and  are 
transmitted  at  boosted  power  level.  Pilots  can  be  used  for  frame  synchronization, 
frequency  synchronization,  time  synchronization,  channel  estimation,  transmission  mode 
identification  and  can  also  be  used  to  compensate  for  hardware  limitations,  such  as  phase 
noise  [9]. 

Data  and  overhead  subcarriers  are  periodically  allocated  in  time  and  frequency  so 
that  channel  variations  can  be  tracked  efficiently  with  minimal  loss  of  data  bandwidth. 
Frequency  allocations  of  data  and  overhead  subcarriers  are  shown  in  Figure  9.  The 
subcarriers  reserved  for  overhead  are  subdivided  into  pilots  (continual  and  scattered),  P2 
symbols  for  defining  various  parameters  used  for  the  transmission  and  frame-closing 
symbols  for  quick  synchronization.  All  other  subcarriers  contain  data.  Of  particular 
interest  are  the  scattered  pilots,  which  cover  a  sufficient  range  of  frequencies  so  that  the 

channel’s  frequency  response  can  be  easily  determined  by  interpolation. 
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The  locations  of  continual  pilots  are  fixed  from  symbol  to  symbol,  and  those  of 
scattered  pilots  change  from  symbol  to  symbol.  Further  details  on  the  distribution  of 
pilots  are  provided  in  [9].  One  of  the  pilot  patterns  used  in  DVB-T2  signal  is  shown  in 
Figure  9. 


P2  symbols 


Normal  symbols 


Frame  closing  symbol 

■  Continual  Pilots  ■P2  Pilots 

■  scattered  Pitots  ■  Frame  Closing  Pilots 

Figure  9.  Pilot  pattern  used  in  DVB-T2. 

7.  Implementation 

DVB-T2  includes  many  configurable  parameters,  as  discussed  previously.  A 
Simulink  model  was  developed  for  the  modulator  and  demodulator  functions  using  the 
configuration  discussed  next.  Each  superframe  consists  of  four  T2  frames,  and  each  T2 
frame  consists  of  68  OFDM  symbols.  In  particular,  in  what  is  called  the  “2K  mode,”  for 
each  OFDM  symbol  there  are  1,512  data  subcarriers,  343  nulls,  176  pilot  subcarriers,  and 
17  P2-subcarriers,  for  2,048  subcarriers.  The  distribution  of  pilot  subcarriers  changes 
cyclically  from  OFDM  Symbol  0  to  OFDM  Symbol  67.  Scattered  and  continual  pilots 
overlap  at  ten  locations  to  make  176  pilots. 

The  GI  length  varies  with  the  FFT  size.  For  example,  for  the  2K  mode,  GI  choices 
are  1/32,  1/16,  and  1/4  [9].  In  the  case  of  a  2048-point  FFT  with  GI  choice  of  1/32, 
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presented  in  this  research,  each  OFDM  symbol  has  a  total  of  2048  +  2048/32  =  2112 
samples.  The  modulation  scheme  for  data  is  QPSK,  and  M-QAM,  with  M  varying  from 
16  to  256.  For  our  simulations,  QPSK  was  chosen  for  data,  while  pilot  subcarriers  are 
transmitted  using  BPSK.  The  block  diagram  of  the  DVB-T2  structure  implemented 
during  simulation  is  shown  in  Figure  10. 

The  proposed  demodulation  scheme  is  based  on  the  detection  and  demodulation 
of  the  continual  and  scattered  pilots  for  channel  estimation  and  the  frame-closing  symbols 
for  fast  frame  synchronization.  This  process  is  an  integral  part  of  our  model-based 
demodulation  using  multiple  antennas  and  Kalman  filtering  and  is  presented  next.  Any 
other  overhead  information  regarding  broadcast  data  demodulation  is  irrelevant  and 
disregarded. 

Underlying  principles  of  the  DVB-T2  standard  were  presented  in  this  chapter  and 
estimation  of  transmission  CIR  characteristics,  using  pilot  signals  present  within  DVB-T2 
signal,  was  also  discussed.  The  approach  proposed  in  this  research  uses  a  two-antenna 
model  and  requires  independence  between  the  two  transmission  channels  to  extract  the 
weak  signal;  therefore,  CIRs  for  both  transmission  channels  are  estimated  from  pilot 
signals  contained  in  the  respective  DVB-T2  signals  received  at  each  antenna.  These  two 
estimated  CIRs  are  then  used  in  a  KF  framework  to  estimate  the  state  of  DVB-T2  signal, 
as  discussed  next  in  Chapter  III. 


Figure  10.  Block  diagram  of  DVB-T2  implementation  in  Simulink. 
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III.  KALMAN  FILTER  APPLICATION  TO  MULTI- ANTENNA 

SYSTEMS 

The  KF  and  its  variants  are  among  the  most  popular  estimation  techniques  and 
have  applications  in  almost  all  fields  of  technology.  This  success  is  due  to  its  simple 
structure,  small  computational  requirements,  recursive  nature,  robustness,  and  optimal 
estimation  for  linear  systems  with  Gaussian  noise.  Estimation  and  different  estimation 
schemes  are  first  discussed  in  this  chapter.  Next,  the  state-space  model  for  a  linear 
dynamic  system  is  described  followed  by  the  application  of  the  KF  to  this  model  to 
estimate  the  state  of  the  system.  Finally,  application  of  the  KF  framework  to  estimate  the 
state  of  DVB-T2  signal  is  discussed. 

A.  ESTIMATION 

Estimation  is  the  process  of  extracting  the  value  of  a  parameter-of-  interest  (POI) 
from  a  random,  noisy,  and  uncertain  data  set  of  measurements.  Noise  can  arise  from 
sensor  measurements,  instrument  inaccuracies,  signal  propagation  channel  or  human 
error  and  causes  randomness  in  POI  as  well.  The  estimator  attempts  to  predict  the  best 
estimate  of  the  unknown  POI  using  these  noisy  measurements.  Filtering  is  the  process  of 
eliminating  noise  (undesirable  signals)  from  noisy  measurements  in  order  to  obtain  the 
best  estimate  of  the  POI.  The  extraction  of  the  best  estimate  requires  some  sort  of 
decision  making  according  to  some  preset  criteria  in  order  to  choose  the  best  from  a  set  of 
alternatives.  The  criteria  can  be  different  for  different  situations,  and  an  optimal  estimator 
yields  the  best  estimate  by  optimizing  certain  criteria.  POI  can  be  a  random  variable, 
which  is  time-invariant,  or  the  state  of  a  dynamic  system,  which  is  a  time-variant 
stochastic  process.  Subsequently,  there  are  two  types  of  estimators  [8]: 

1.  Parameter  estimators  (static  estimation  problem);  for  example,  Weiner- 
Hopf  filter 

2.  State  estimators  (dynamic  estimation  problem);  for  example,  KF 
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B.  ESTIMATION  TECHNIQUES 

The  following  estimation  techniques  can  be  employed  if  the  POI  is  a  random 
variable: 

•  maximum  a  posteriori  (MAP)  estimator 

•  minimum  mean-square  error  (MMSE)  estimator 

•  linear  minimum  mean-square  error  (LMMSE)  estimator 

It  turns  out  that  the  KF  is  the  MMSE  state  estimator  if  the  initial  state  and  noise 
entering  into  the  system  are  Gaussian.  In  the  non-Gaussian  case,  the  KF  is  still  a  LMMSE 
state  estimator  [10].  To  get  a  good  understanding  of  the  KF  problem  and  its  solution, 
some  underlying  concepts  related  to  KF  are  discussed  next. 

C.  THE  KALMAN  FILTER 

Let  x(^)=[Xi  (^),X2  (A:),....,x„  (A:)]^  be  the  state  of  a  discrete-time,  linear  system 

and  u(k)=[Uj  (k),U2  (k),....,u^  (k)]^  be  the  deterministic  input  to  that  system  at  time  k. 

Then  the  system  state  at  time  A:  + 1  is  given  by  the  state  equation 

x(A:  +  l)  =  F(A:)x(A:)  +  G(A:)u(A:)  +  y(A:)  A:  =  0,1, . ,  (3.1) 

where  F(A:)  is  the  x  state  transition  matrix,  G(A:)  is  the  x  input  gain  matrix, 
and  y(A:) is  a  zero  mean  white  Gaussian  process  noise,  with  covariance  matrix 

^{v(A:)y(A:)'}  =  Q(A:). 

If  z(A:)=[z,  (A:),Z2  (A:),....,z^  (A:)]^  is  the  measurement  vector,  then  the 
measurement  or  output  equation  is  given  by 

z(A:)  =  H(A:)x(A:) -I- ra(A:)  A:  =  1,2,....,  (3-2) 

where,  H(A:)  is  the  x  measurement  matrix  and  ro(A:)  is  x  1  zero  mean  white 
Gaussian  measurement  noise  with  covariance  matrix  R(A:)  =  E  {ro(A:)oj(A:)'} . 

The  underlying  assumptions  of  the  previous  linear  dynamic  model  are: 
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1. 


Sampling  interval  is  constant 


2.  Matrices  F(A:),  G(A:)  and  H(A:)  and  deterministic  input  u(k)  are  known  for 
all  values  of  k, 

3.  Yhe  noise  terms  are  zero  mean,  white,  and  uncorrelated,  which  leads  to 
E  {y(A:)®(/)'}  =  0  for  all  values  of  k,j. 

The  notation  k.{k\k)  represents  the  estimate  of  Mk)  at  time  k  based  on  all 
measurements  up  to  and  including  z{k) .  Similarly,  x{k  +  \\k)  represents  the  estimate  of 
x{k  + 1)  at  time  k  + 1  based  on  all  measurements  up  to  and  including  z{k) .  This  notation 
is  used  extensively  in  the  follow-up  discussion.  The  symbols  P^^{k\k^  and/^ (A:-l-l|A:) 
represent  the  covariance  matrices  of  state  estimation  errors  at  time  k  and  k  +  \, 
respectively.  The  KF  estimates  the  state  of  the  system  when  a  new  measurement  is 
obtained.  Updating  the  state  estimate  when  a  new  measurement  is  taken  can  be  split  into 
two  steps,  prediction  and  correction,  which  are  presented  next. 

1.  Prediction 

State  estimate,  state  covariance,  and  measurement,  respectively,  are  predicted  to 
the  next  measurement  time  using  the  following  equations: 

k{k  +  \\k)  =  Kk)k{k\k)  +  G(k)u(A:),  (3.3) 

F{k  +  l\k)  =  F(k)Pik\k)F{ky+Q{k),  (3 .4) 

and 

z{k  +  l\k)  =  U(k  +  l)x{k  +  l\k).  (3.5) 

State  prediction  error  and  measurement  prediction  error  are  given  by 

x(A:  -I-  \\k)  =  x(A:+l)  -  x(A:  -I-  \\k)  =  F{k)x{k\k)  +  yik)  (3-6) 

and 

z(k  +  1|A:)  =  z{k  + 1)  -  z{k  +  \\k) 

=  H(A:  +  l)x(A:  +  l|A:)  +  w(A:  +  l),  (3.7) 

respectively,  where  l{k  +  \\k)\?.  called  innovation  or  measurement  residual.  The 
measurement  prediction  (innovation)  covariance  matrix  S  is  computed  as 
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S(k  + 1)  =  n{k  + 1)  P(A:  +  \\k)n{k  + 1)'  +  R(A:  + 1), 


(3.8) 


and  the  filter  gain  is 

W(it  + 1)  =  ?(k  +  \\k)li(k  +  \yS(k  + 1)"‘ .  (3.9) 

2.  Correction 

The  updated  state  estimate  and  state  eovariance  at  time  step  k  +  l  are  given  by  the 
state  update  equation  and  covariance  update  equation,  respectively,  as  [10]: 

x(A:  +  l|A:  +  l)  =  x(A:  +  l|A:)  +  W(A:  +  l)z(A:  +  l|A:)  (3.10) 

and 

P(A:  + 1|  A:  + 1)  =  P(A:  + 1|  A:)  -  ?{k  +  \\k)ll{k  +  \yS{k  + 1)"’  H{k  +  1)P(A:  +  l|yi:) 

=  P(k  + 1|  A:)  -  W(k  +  \)S(k  +  l)W(k  +  1)'.  (3.11) 

The  flow  chart  of  one  KF  cycle  is  presented  in  Figure  1 1 .  The  state  estimation  cycle 
consists  of  state  and  measurement  prediction  (time  update)  and  state  update 
(measurement  update).  The  covariance  calculations  are  independent  of  the  state  and 
measurements  and  can  be  performed  offline. 

The  computational  requirements  of  the  KF  are  approximately  proportional  to 
n^,n  =  max(n^,n^) .  In  most  cases,  a  steady-state  implementation  in  which  the  KF 

parameters  (state  error  covariance  matrix  and  Kalman  gain)  are  kept  constant  and  pre¬ 
computed  yields  satisfactory  results  with  considerable  savings  in  computations. 

3.  Riccati  Equation 

The  expression  for  the  state  covariance  matrix  P,  given  in  Equation  (3.4),  can  be 
written  as 

Pik  +  1|A:)  =  F(A:)  {P(k\k  - 1)  -  Pik\k  -  1)H(A:)'  S(A:)H(A:)P(A:|A:  - 1)}  x 

F(A:)'+Q(A:).  (3.12) 

Equation  (3.12)  is  known  as  the  algebraic  matrix  Riccati  equation,  and  its  solution 
for  a  LTI  system  converges  to  a  finite  steady-state  (finite)  covariance  matrix  if  the  pair 
{F,H}  is  completely  observable;  i.e.,  the  state  can  be  fully  determined  from  the  input  and 
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output  measurements.  The  steady-state  covariance  matrix  is  the  solution  of  the  algebraic 
matrix  Riccati  equation,  or  simply  the  algebraic  Riccati  equation,  and  yields  the  steady- 
state  gain  for  the  KF,  which  is  given  as 

K  =  FPH'iHPH' +  Ry\  (3.13) 


Evolution 
of  the  system 
(true  state) 


Known  input 

(control  or  Estimation 

sensor  of  the  state 

motion) 


State  covariance 
computation 


Figure  11.  One  cycle  of  KF  state  estimation,  from  [8]. 


D.  KALMAN  FILER  ESTIMATION  OE  MULTI-ANTENNA  SYSTEMS 

Let  x{n\  be  the  transmitted  DVB-T2  signal  and  .s[«]  be  the  weak  signal  that  is  to 
be  recovered.  The  bandwidth  of  5[«]  is  assumed  to  be  much  smaller  than  that  of  x[n\ , 
with  a  low  data  rate  and  longer  symbol  duration  than  that  of  x{n\ .  Such  characteristics 
ensure  that  the  time  spread  of  the  two  channels  is  negligible  as  compared  to  the  symbol 
duration  of  .s[n],  and  multipath  can  be  neglected  with  the  two  antennas  receiving  the 

same  copy  of  5'[n] . 
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The  signals  received  at  the  two  antennas  are  given  by 

y\n\  =  h\n\*x{n\  + s{n\  +  w^{n\  (3.14) 

and 

y2[n\  =  h2{n\*x{n\  + s{n\  +  W2{n\,  (3.15) 

where  h\n\  and  h2\n\  are  CIRs  that  are  time-invariant  and  independent  of  each  other; 
Wj[n]  and  ^^[n]  are  zero  mean  additive  white  Gaussian  noise  at  the  two  antennas  with 

covariance  matrices  R  =  £'|wj[m]wj[ny| ,  and  0  =  £'|n’2[m]w2[n]^| .  We  also  assume  the 

noise  sequences  to  be  uncorrelated  with  each  other,  which  leads  to  £'|wj[m]w2[n]^|  =  0 
for  all  values  of  m  and  n  . 

Equations  (3.14)  and  (3.15)  can  also  be  written,  respectively,  as 

fi  [^]  =  +  h^\n\x{n  - 1]  -i- +  h^j^{n\x{n  -  A^]  -i-  ^[n]  -i-  w\n\  (3.16) 

and 

y2\n\  -  h2f\n^x{n\  +  h2^\n\x\n-\\  + +  h2i^{n\x\n-N'\  +  s\n\  +  w^n\.  (3.17) 


Rewriting  Equation  (3.16),  we  get 


An\  =  — —  {Ti  [«]-/?!  1  [n\x[n  - 1]  - -  h^2,[n\x[n  -N]-  s[n\  -w\_n\}.  (3.1 8) 

KM 


Next,  Equations  (3.18)  and  (3.17)  can  be  written  in  matrix  form,  respectively,  as 

x[«]  ~KlK  ~KlK  ~Kn-\)IK  ~Kn/K  x[«— 1] 

x[«-l]  1  0  •••  0  0  x[«-2]  ^ 

x[n-(A^-l)]_  0  O'"  1  0  _x[n-A^]_ 

VK  ~yK 

;  T,W+  ;  [s[n]  +  w,[n]]  (3.19) 

0  J  L  0  _ 

and 
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y2W  =  [^2i  Ki  ••• 


^[s[n\  +  w^[n\]. 


(3.20) 


x[n  - 1] 
x[n-N} 


where  has  been  forced  to  zero;  i.e.,  h^in]  is  delayed  by  one  time  step. 


Defining  x[n]  =  {x[«-l],x[n-2], . ,x[«-A^]}^,  we  can  write  Equations  (3.19) 

and  (3.20),  respectively,  in  state  space  form  as 

x[n  +  l]  =  ^x[n]  +  5yj[n]  +  G{5[n]  + Wj[n]}  (3.21) 

and 

y2\n\  =  C^x[n]  +  {^[n]  +  >V2[n]} ,  (3.22) 


where  the  matrices  A,  B,  G  and  C  are  defined,  respectively,  as 

f— /?ll//?10  “^12/^10  ■■■  ~K(N-Y)IKo  ~\n/\o 


A  = 

1 

0 

0 

0 

0 

1 

'-VK 

^20 

B  = 

0 

,  G  = 

0 

,  C  = 

^21 

_  0 

0 

J^IN  _ 

(3.23) 


(3.24) 


Equations  (3.21)  and  (3.22)  yield  the  state-space  model  which  was  described  by 
Equations  (3.1)  and  (3.2).  The  recursive  estimate  of  the  state  x[n]  can  be  computed  by 

the  KF,  where  the  received  sequence  y\n\  is  the  known  KF  input  and  the  received 

sequence  y2[^]  the  measurement. 


Since  neither  CIRs  change  over  time,  the  matrices  A,  B,  C  and  G  are  constant. 
Writing  Equations  (3.10)  and  (3.11)  according  to  the  state-space  model  of  Equations 
(3.21)  and  (3.22),  we  get 

x[n  + 1]  =  Ax{n\  +  By,  {n\  +  i:[n]  {y^  [n]  -  Cx[n]} .  (3 .25) 
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At  this  point,  Equation  (3.23)  can  be  rewritten  in  the  following  form 
x[n  + 1]  -  -  KC^  x[n]  =  By\n\  +  Ky^iril, 


(3.26) 


where  K  is  the  KF  gain  in  steady-state  and  is  given  by  Equation  (3.13).  Rewriting 
Equation  (3.24)  in  the  z-domain  leads  to 

Xiz)  =  \{zI-{K-AC)y'  b\  Y^{z)  +  \{zI-{K-AC)y'  k\  Y^iz),  (3.27) 

where  X(z)  =  [z"'2r(z),z""X(z),---,z"^2r(z)]^  is  the  z-transform  of  the  sequence 
x[n]  =  {x[n  — l],x[n  —  2], . ,x[n  -  . 

It  can  be  shown  that  for  x[«]  =  {x[n-l],x[«-2], . ,x[n-A]}^, 

x[n-l]  =  [l  0  0---]x[n],  x[n-2]  =  [0  1  0---]x[n]  and  so  on;  therefore,  estimates  of 
states  at  time  n  - 1  and  time  n-2  are  given,  respectively,  by 

z-’X(z)  =  [l  0  0---][{z/-(X-^C)}”‘r  Y^(z)  + 

[1  0  0---'^{zI-iK-AC)y' K  Y^Cz)  (3.28) 

and 

z-"X(z)  =  [0  1  0---'^{zI-{K-AC)y' B  l^(z)  + 

[0  1  0---][{z/-(X-^C)}"‘i:  Y^{z).  (3.29) 

Similarly,  estimates  of  states  at  other  time  steps  can  be  found  using  the  logic 
described  before. 

From  Equations  (3.26)  and  (3.27),  it  can  be  concluded  that  two  transfer  functions 
are  needed  to  produce  the  estimate  of  transmitted  data  at  different  discrete  times  when 
applied  to  the  received  sequences  y\n\  and  y2[«];  thus,  estimating  the  state  x[n-l] 
requires  two  transfer  functions 

Gi(z)  =  [l  0  0---][{z/-(X-^C)}'’r  (3.30) 

and 


26 


(3.31) 

(3.32) 

(3.33) 

(3.34) 

(3.35) 


Equation  (3.33)  can  also  be  written  as 

x{n-\]  =  J\\]y\n]  +  J\l]y\n-\\  +  --- 

2 [1]3^2 M  +  *^2 [2] 3^2 1]  - 

-v[l]x[n-2]-v[2]i[n-3] -  (3.36) 

In  matrix  form,  we  can  write  Equation  (3.34)  as 

yx\n\  1  \y2{n\ 

x[n-l]  =  [J,[l]  J,[2]-]  y,[n-l]  +[/,[!]  J,[2]-]  y,[n-l]  + 

x[n  -  2] 

[-v[l]  -v[2]--]  x[«-3]  .  (3.37) 

Equation  (3.35)  is  implemented  to  estimate  the  state  of  the  transmitted  DVB-T2 
signal.  The  Kalman  gain  parameter  is  calculated  using  the  MATLAB  function  ‘dlqe.m’, 
while  coefficients  of  the  two  transfer  functions  are  calculated  using  the  MATLAB 
function  ‘ss2tf  m’  function.  These  transfer  functions  are  implemented  using  two  discrete 
time  filters,  and  received  sequences  are  then  passed  through  these  two  filters.  Next,  the 
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outputs  are  combined  to  provide  the  DVB-T2  signal  estimate.  This  process  is  illustrated 
in  Figure  12.  Note  that  Kalman  gain  and  coefficients  of  the  transfer  functions  are 
calculated  offline  at  the  beginning  of  the  simulation  and  then  updated  only  periodically. 


data. 

In  this  chapter,  we  discussed  the  KF  based  approach  designed  to  use  the  two 
estimated  CIRs  and  signals  received  at  the  antenna  to  estimate  the  state  of  the  strong 
interferer  DVB-T2  signal.  Specifically,  we  showed  that  we  can  obtain  the  estimate  of  the 
signal  received  at  the  second  antenna  by  using  the  estimated  DVB-T2  signal  and 
estimated  CIR  of  the  second  channel.  Subtraction  of  estimated  signal  from  the  original 
signal  received  at  second  antenna  suppresses  the  interference,  which  allows  the  recovery 
of  the  weak  signal.  This  approach  was  implemented  in  Simulink,  and  simulation  results 
are  discussed  in  Chapter  IV. 
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IV.  SIMULATION  AND  RESULTS 


Theoretical  background  to  recover  a  weak  signal  in  the  presence  of  a  strong  co¬ 
channel  interferer  was  presented  in  preceding  chapters.  Simulink  models  were  developed 
to  evaluate  the  robustness  of  the  proposed  approach  to  noise  distortions.  In  this  chapter, 
we  first  present  the  scheme  implemented.  Next,  we  introduce  concepts  of  SNR  and  SIR. 
Finally,  we  investigate  the  approach  robustness  to  SNR,  SIR,  and  estimated  CIR 
accuracy. 

As  mentioned  earlier,  the  received  signal  at  each  antenna  consists  of  a  strong 
DVB-T2  signal  (after  passing  through  a  channel),  a  weak  signal-of-interest  and  additive 
white  Gaussian  noise.  In  this  study,  the  weak  signal  5'  is  a  narrow-band  signal  modulated 
to  a  carrier  and  is  embedded  in  interference  and  noise  as  shown  in  Figure  13. 
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Figure  13.  Frequency  domain  plot  of  interferer  (>>2) ,  weak  signal  (s) 
and  noise  (w)  for  SNR=20  dB  and  SIR=-25  dB. 


It  has  been  shown  theoretically  that  the  state  of  a  signal  can  be  estimated  using  the 
KF  framework  described  in  Chapter  III.  The  said  framework  is  used  to  estimate  the 
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transmitted  DVB-T2  sequence,  which  requires  knowledge  of  CIRs.  Pilots  present  in  the 
DVB-T2  signal  are  used  to  estimate  both  CIRs.  Convolution  of  the  estimated  DVB-T2 
sequence  with  one  of  the  estimated  CIR  provides  us  with  an  estimate  of  the  signal 
received  at  the  antenna  corresponding  to  that  channel.  The  subtraction  of  the  estimated 
TV  broadcast  signal  from  the  overall  signal  received  at  the  antenna  is  designed  to  remove 
the  contribution  of  the  strong  TV  broadcast  signal  and  preserve  the  weak  signal  and 
noise,  from  which  the  weak  signal  can  be  extracted.  This  scheme  is  presented  in 
Figure  14. 


Figure  14.  Block  diagram  for  recovering  weak  signal. 


As  evident  from  the  above  scheme,  the  strong  interferer  contribution  can  be 
minimized,  as  shown  in  Figure  15,  provided  accurate  estimates  of  x ,  ^  and  are 

obtained.  The  resulting  signal  is  a  combination  of  the  weak  signal  and  noise  which  then 
becomes  a  standard  demodulation  problem,  and  the  weak  signal  can  be  recovered  easily, 
as  illustrated  in  Figure  16. 
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Figure  15.  Suppression  of  interferenee  and  reeovery  of  weak  signal  for 

SNR=20  dB  and  SlR=-25  dB. 


Frequency  (MHz) 

Figure  16.  Reeovered  weak  signal  (s)  after  demodulation  for 


SNR=20  dB  and  S1R=  -  25  dB. 
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A.  SIMULINK  IMPLEMENTATION 

The  scheme  shown  in  Figure  14  was  implemented  in  Simulink,  using  x  and  5  as 
zero  mean  QPSK  signals.  The  signal  x  is  implemented  at  baseband,  while  s  is  band 
limited  by  a  square  root-raised  cosine  fdter  and  modulated  at  a  carrier  frequency  equal  to 
914.286  kHz.  Since  both  x  and  s  are  complex  signals,  zero-mean  independent,  complex 
white  Gaussian  noise  and  are  added  at  both  antennas. 

For  simplicity,  no  coding  was  applied  and  perfect  synchronization  both  in  time 
and  frequency  was  assumed.  The  results  presented  further  in  the  research  were  obtained 
by  running  the  simulation  for  10^  symbols  of  the  weak  signal.  Performances  are  evaluated 
in  terms  of  bit  error  ratio  (BER)  of  both  DVB-T2  and  weak  signals  at  different  SNRs  and 
SIRS. 

1.  Signal-to-Noise  Ratio  and  Signal-to-Interference  Ratio 

Recall  that  two  sources  of  disturbances  are  present:  the  TV  broadcast  interference 
and  the  AWGN.  Thus,  two  power  level  ratio  quantities  relative  to  the  signal-of-interest 
are  the  signal-to-interference  ratio  and  the  signal-to-noise  ratio,  respectively,  defined  as 

SIR  =  %  and  SIR,s  =  20  logi^  —  ,  (4.1) 

where  errand  a  represent  variance  and  standard  deviation,  respectively,  and  x  is 
normalized  to  1,  which  means  cr^  =  cr^  =  1 .  Thus,  for  a  given  SIR,  cr^  is  given  by 

SIR 

0-^=10  20.  (4.2) 

Since  s  is  band  limited  with  bandwidth  equal  to  91.2  kHz  and  the  noise  is  white, 
the  term  SNR  can  be  defined  as 

5VR,,=201og,j^ 

f  \ 

=  201ogjo  +201ogjo 

V  ) 
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where  is  the  two-sided  noise  power  spectral  density,  and  and  represent  the 
bandwidths  of  weak  signal  -y  and  white  noise  w,  respectively.  For  a  given  SNR  and  <t^  , 
cr^  can  be  obtained  from  Equation  4.2  as 

SNR-20log^„{B„lB,)/ 

C7^  =  CT,10  /20.  (4  4) 

Note  that  both  CIRs  are  normalized,  meaning  the  power  of  the  signals  passing 
through  the  channels  remains  the  same.  Since  =  I,  the  signals  Pj  and  yj 

received  at  the  antennas  also  have  variances  =  cr^  =  1 .  Therefore,  the  preceding 

yi  yi  A 

definitions  of  SNR  and  SIR  also  hold  for  signals  Pj  and  y2 .  For  simulation  purposes,  we 

first  defined  SIR  and  SNR  and  set  standard  deviations  for  s  and  w  to  achieve  the  desired 
SIR  and  SNR. 

2.  Estimated  CIR 

Recall  that  DVB-T2  pilots  are  used  to  estimate  the  CIRs  for  equalization  in  DVB- 
T2  receiver  and  to  estimate  the  state,  as  shown  earlier  in  Figure  14.  Simulation  results 
showed  that  BER  of  10”^  and  below  were  achieved  for  the  estimated  DVB-T2  signal 
using  estimated  CIRs;  however,  the  same  type  of  performance  was  not  observed  for  the 
weak  signal  BER  when  the  estimated  CIR  parameters  were  used  in  the  KF  step  for  the 
DVB-T2  signal  state  estimation.  In  this  case,  the  best  BER  achieved  for  the  weak  signal 
was  equal  to  0.15;  hence,  results  show  that  estimating  the  CIR  characteristics  using  the 
pilots’  information  is  sufficient  for  DVB-T2  signal  demodulation  but  is  not  for  the  weak 
signal  demodulation  process.  Therefore,  original  CIRs  coefficients  and  variants  with  5%, 
10%,  and  15%  random  errors  in  coefficient  values  were  used  as  estimated  CIRs  in  the 
simulations  to  evaluate  the  sensitivity  of  the  overall  proposed  scheme  to  degradations  in 
estimated  CIR  characteristics. 

Results  showed  that  the  estimation  process  breaks  down  when  CIR  estimates  have 
a  20%  error.  At  that  point,  both  estimated  CIRs  carry  the  same  information,  and  the 
independence  condition  of  the  two  channels  needed  for  the  scheme  to  work  breaks  down. 
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In  turns,  the  KF  implementation  becomes  ill-conditioned  and  is  unable  to  predict  the  state 
of  the  DVB-T2  signal. 

B.  SIMULATION  RESULTS 

The  simulation  scheme  mentioned  above  was  implemented  in  Simulink.  The 
results  obtained  are  discussed  below. 

1.  Signal  Attenuation 

As  noted  above,  recovering  the  weak  signal  is  possible  only  when  the 
transmission  channels  associated  with  the  two-antenna  receiver  are  independent,  in  the 
sense  that  their  associated  transfer  functions  have  different  root  locations,  ensuring  the 
combined  state  space  system  is  fully  observable  and  controllable. 

We  investigated  the  effectiveness  of  this  system  by  observing  the  attenuations  of 
the  two  received  broadcast  signals  Pj  and  >^2 »  the  signal-of-interest  s  and,  finally,  the 

noise  w .  Since  the  parameters  at  the  receiver  are  adjusted  on  the  basis  of  the  estimated 
multipath  of  the  TV  broadcast  signal,  sensitivity  of  these  parameters  to  errors  in  the 
channel  estimates  is  also  inspected. 

Recall  that  the  broadcast  TV  signals  Pj  and  >^2 »  received  through  the  multipath 

channels,  need  to  be  highly  attenuated  with  respect  to  the  signal-of-interest  s  while  the 
noise  w  is  not  amplified  by  the  receiver  for  the  proposed  scheme  to  be  successful.  Thus, 
we  tested  the  receiver  ability  to  separately  attenuate  signals  Pj ,  y2,  s,  and  w  by 
broadcasting  the  following  signal  combinations: 

•  DVB-T2  signals  Pj  and  >^2  only 

•  noise  w  only 

•  signal-  of-interest  s  only 

In  each  case,  we  computed  the  power  attenuation  (or  gain)  level  in  dB  between 
each  test  signal  (Pjand  y2,s  ,w)  as  input  and  the  resulting  output  (>’2  =  y2-K)-  The 
attenuation  parameter  is  defined  as 
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(4.5) 


Attn^{dB)  =  lOlogjo 

where  cr^  represents  the  variance  of  the  signal  whose  attenuation  level  is  to  be  observed. 
a.  Interference  Attenuation 

First  we  apply  the  broadcast  TV  signals  and  y2  as  inputs,  without  noise  w  or 
signal-of-interest  s  present.  Attenuation  levels  for  signals  y^  and  yj  obtained  with 

actual  CIRs  and  with  5%,  10  %,  15%  error  levels  for  estimated  CIRs  are  plotted  in  Figure 
17  and  Figure  18,  respectively.  Results  show  that  when  actual  CIRs  are  used  as  estimated 
CIRs,  signals  y^  and  y2  are  completely  suppressed.  The  level  of  attenuation  obtained  for 

signals  y^  and  yj  increases  with  the  level  of  accuracy  in  the  estimated  CIR  parameters.  It 
can  also  be  observed  that  the  attenuation  level  is  insensitive  to  noise  power  changes. 
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Figure  17.  Attenuation  of  the  interference  y^  for  different  noise  power 
levels  and  error  percentages  in  the  estimated  CIR  coefficients. 


It  is  worth  noting  that  the  KF  becomes  deterministic  when  exact  CIR  parameters 
are  known  and  no  other  disturbance  (noise  w  or  signal-of-interest  5)  is  present; 
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therefore,  the  innovation  of  the  KF  is  equal  to  zero  in  steady-state.  Since  this  case 
represents  a  deterministic  observer,  the  state  estimation  error  in  steady-state  is  equal  to 
zero,  which  leads  to  infinite  attenuation. 


A — Actual  CIRs 
— 5%  error 
B — 10%  error 
— 15%  error 


Figure  18.  Attenuation  of  the  interference  y2  for  different  noise  power 
levels  and  error  percentages  in  the  estimated  CIR  coefficients. 


Any  other  CIR  parameter  error  yields  nonzero  attenuation  for  signals  and  Pj  • 

Due  to  the  optimality  of  the  KF,  the  sensitivity  to  errors  in  the  estimated  CIR  parameters 
is  minimized,  which  is  indicated  by  a  high  attenuation  (about  50  dB)  for  the  range  of 
noise  variance  levels  considered  in  this  implementation. 

b.  Signal  of  Interest  Attenuation 

In  the  second  case,  we  apply  the  signal-of-interest  s  as  input  without  broadcast 
TV  signals  Pj  and  ^2  or  noise  w  present.  Attenuation  levels  for  s  obtained  with  actual 
CIRs  and  with  5%,  10  %,  and  15%  error  levels  for  estimated  CIRs  are  plotted  in  Figure 
19.  It  can  be  seen  from  Figure  19  that  s  is  attenuated  by  approximately  9.5  to  10  dB,  and 
the  level  of  attenuation  observed  for  s  is  independent  of  noise  power  level  variations  in 
the  KF. 
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Figure  19.  Attenuation  of  the  signal  s  for  different  noise  power  levels 
and  error  percentages  in  the  estimated  CIR  coefficients. 

c.  Noise  Attenuation 

In  the  third  case,  we  apply  the  noise  w  as  input  without  broadcast  TV  signals 

and  y2  or  signal-of-interest  s  present.  Attenuation  levels  for  w  obtained  with  actual 

CIRs  and  with  5%,  10  %,  and  15%  error  levels  for  estimated  CIRs  are  plotted  in  Figure 
20.  Results  show  the  noise  w  attenuation  is  about  10  to  10.5  dB.  Again,  the  attenuation 
level  does  not  change  significantly  with  noise  power  levels  variations. 

It  has  been  shown  that  y^  and  y2  suffer  far  more  attenuation  than  that 
experienced  by  5 ,  and  ;  hence,  3)2  consists  of  predominantly  the  weak  signal  s 
and  the  noise  w  and  only  includes  a  negligible  contribution  from  the  signal  y2 . 

A  basic  requirement  of  the  system  considered  is  that  the  signal-of-interest  has  low 
power  so  that  it  does  not  interfere  with  TV  broadcasting  demodulation.  Having 
established  that  the  scheme  presented  in  the  research  suppresses  the  interference  in  all 
noise  conditions,  we  now  evaluate  its  performance  by  calculating  BER  for  the  DVB-T2 


37 


and  the  weak  signal  at  different  SIRs  and  SNRs.  Our  goal  is  to  identify  some  SIR  and 
SNR  where  BER  is  less  than  10“^  in  order  to  correctly  demodulate  both  the  interference 
and  the  weak  signal.  For  this  reason,  BER  is  plotted  against  SNR  at  different  SIRs  for  the 
DVB-T2  signals,  received  from  two  different  channels,  and  the  weak  signal. 


Figure  20.  Attenuation  of  the  noise  w  for  different  noise  power  levels 
and  error  percentages  in  the  estimated  CIR  coefficients. 

2.  DVB-T2  BER  Performance  Analysis 

It  is  evident,  as  expected,  from  Figure  21  and  Figure  22  that  the  BER  obtained  for 
the  DVB-T2  signal  decreases  with  smaller  SIR. 

For  the  signal  transmitted  through  channel  I,  a  BER  below  10“^  was  obtained 
with  SNR=10  dB  and  SIR=  -  25  dB.  However,  for  the  signal  transmitted  through  channel 
2,  a  BER  of  less  than  10“^  was  achieved  with  SNR=25  dB  and  SIR=-25  dB.  Hence, 
these  results  show  that  the  BER  for  the  DVB-T2  signal  depends  on  the  noise  level  present 
in  the  channel.  Some  type  of  signal  processing  should  be  considered  to  select  the  signal 
from  the  channel  which  has  the  lower  BER. 
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Figure  21.  DVB-T2  Channel  1  signal  BER  at  different  SIR. 


Figure  22.  DVB-T2  Channel  2  signal  BER  at  different  SIR. 
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3.  Weak  Signal  BER  Performance  Analysis 

Due  to  the  high  attenuation  of  TV  interference  in  the  proposed  system,  the 
contribution  of  the  TV  component  in  the  resulting  signal  (^2  =^2  “A)  is  smaller  than 

that  of  the  signal-of-interest  and  noise.  Consequently,  results  show  that  the  BER  for  the 
weak  signal  is  not  affected  by  variations  in  the  TV  signal  power  and,  therefore,  BER 
remains  constant  for  all  SIR  considered. 

It  was  shown  earlier  that  the  attenuation  level  of  the  broadcast  TV  interference 
increases  as  estimated  CIRs  coefficients  accuracy  improves;  however,  this  had  no  bearing 
on  the  BER  of  the  weak  signal.  Results  indicated  that  with  0  to  15%  errors  in  the 
estimated  CIR  coefficients,  the  BER  obtained  for  the  weak  signal  remained  constant. 
When  estimated  CIR  coefficients  have  20%  or  more  error,  the  KF  estimation  process 
breaks  down,  and  the  weak  signal  cannot  be  recovered. 


Figure  23.  Weak  signal  BER  with  SIR=  -  5  dB  and  15%  error  in 

estimated  CIR  coefficients. 

Simulations  showed  that  the  BER  obtained  for  the  weak  signal  remain  constant 


irrespective  of  the  SIR  or  when  estimated  CIRs  have  0  to  15  percent  errors  in  their 

40 


coefficients.  Therefore,  only  BER  obtained  for  the  weak  signal  with  SIR=  -  5  dB  and 
15%  error  level  in  estimated  CIRs  coefficients  are  presented  in  Figure  23.  It  is  observed 
from  Figure  23  that  a  minimum  SNR  equal  to  7dB  is  required  to  ensure  a  maximum  BER 
of  10’^ .  Similar  results  were  obtained  at  different  SIRs  and  with  different  error 
percentages  in  the  estimated  CIR  coefficients. 

Overall  results  indicate  that  the  proposed  scheme  can  be  used  to  recover  the  weak 
signal-of-interest  from  the  co-channel  interference  signal  without  distorting  the 
interference  signal  provided  estimated  CIRs  have  15  percent  or  less  errors  in  their 
coefficients  and  the  SIR  is  below  -25  dB .  Results  are  summarized  in  Chapter  V,  and  our 
conclusions  and  recommendations  for  future  work  are  presented. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  purpose  of  this  thesis  was  to  recover  a  weak  signal  in  the  presence  of  a  strong 
interference  signal  and  noise.  The  interference  signal  was  assumed  to  be  affected  by 
multipath  and  the  two  paths  present  to  be  independent  of  each  other.  The  interference 
component  was  modeled  as  a  DVB-T2  signal  at  baseband,  while  the  weak  signal  was 
selected  as  a  narrowband  carrier  modulated  QPSK  signal  not  affected  by  multipath. 

A  two-antenna  model  using  KF  was  presented  to  estimate  the  state  of  the 
interference  signal  using  channel  diversity  at  the  receiver.  The  known  overhead  in  the  TV 
broadcast  interference  signal  was  used  to  estimate  the  transmission  Channel  Impulse 
Responses  (CIRs)  that  are  used  to  equalize  the  TV  broadcast  signal  with  the  resulting 
BER  below  10“^  for  the  DVB-T2  signal.  However,  a  BER  of  10’^  or  below  for  the 
estimated  weak  signal  could  not  be  obtained  as  the  channel  estimation  scheme  based 
using  only  the  DVB-T2  pilots  information  was  not  accurate  enough.  Therefore,  the 
research  conducted  in  this  study  assumes  that  the  transmission  channel  estimation  can  be 
done  accurately  and  investigates  the  robustness  of  the  proposed  scheme  to  estimated  CIR 
coefficients  level  of  accuracy. 

The  proposed  KF  demodulation  is  based  on  the  CIRs,  and  the  approach  is  shown 
to  be  effective  in  the  presence  of  some  level  (up  to  15%)  of  errors  in  the  CIR  estimates. 
Results  also  show  that  the  approach  breaks  down  when  CIR  estimates  have  20%  or  more 
error.  At  that  point,  both  estimated  CIRs  carry  the  same  information,  and  the  requisite 
independence  condition  between  the  two  channels  no  longer  holds  true.  The  KF 
implementation  becomes  ill-conditioned  and  is  unable  to  predict  the  state  of  the  DVB-T2 
signal. 

Results  also  show  that  significant  TV  broadcast  interference  attenuation  is 
obtained  when  using  actual  or  estimated  CIRs  in  the  KF  operation.  In  such  cases,  the 
signal  obtained  afterwards  is  mostly  composed  of  the  weak  signal  and  noise,  out  of  which 
the  weak  signal  can  be  easily  recovered.  Results  obtained  indicate  that  interference  BER 
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improves  with  decreasing  SIR.  In  addition,  it  has  also  been  observed  that  the  weak  signal 
BER  is  independent  of  the  SIR  and  is  a  function  of  SNR  alone. 

Follow-on  research  should  include  a  more  robust  approach  for  estimating  the 
channel  coefficients  so  that  it  can  be  seamlessly  integrated  into  the  existing  model. 


44 


APPENDIX.  MATLAB  CODE 


Matlab  code  used  in  the  research  along  with  associated  functions  is  given  below. 
1.  Main  File 

clc 

clear  all 
close  all 
%%  Load  Model 

model  =  ' Kalman0ne_hlh2_Actual_Latest_SNR ' ; 
load_system (model ) 

stop_weak  =  le3; 
stop_int  =  inf; 
rho  =  0/100; 
iter  =  100; 

%%  Transmission  channel  parameters 

hl=[0.7  +  0.3i,  2. 3-. 11,  -0.9  +  0.41,  0.21]; 
h2=[0.3+0.9i,  0.1-0.41,  -1.31,  2.2+.02i]; 

hi  =  [hi , zeros ( 1 , 10-length (hi ))] ; 
h2  =  [h2 , zeros ( 1 , 10-length (h2 ))] ; 

hi  =  hi  /  sqrt  (hl’^hl  '  )  ; 
h2  =  h2  /  sqrt  (h2’'h2  '  )  ; 

Up  =  100;  %Up  sampling  &  down  sampling  rate  for  weak 

signal 

%  Frequencies 
Fs=  (64/7)  ’^leOG; 
rate  for  DVB-T2  signal 
M=4; 

M_weak  =  4; 
ph_weak  =  pi/4; 

Fb=log2  (M)  (1512/  (2048  +  64)  )  ’^Fs; 

%  Weak  Comms  Signal 

Fc=Fs/10;  %  carrier  frequency  for  weak  signal 

Fdw=Fc/10;  %  data  rate  for  weak  signal 

%%  DVB-T  Settings 

%  DVB-T  settings : 

%  Kdata  =  data  subcarriers  (68x1512) 


%  sampling  and  symbol 

%  M-QAM 
%  M-PSK 

%  phase-offset 

%  bit  rate  for  DVB-T2  signal 
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%  Kpilots  =  pilots  subcarriers  (68x176) 
%  Ktps  =  TPS  subcarriers  (68x17) 

%  1512+176+17=1705 
%  nulls=2048-1705=343 


[Kmod,  Kdem,  Kdata,  Xoh,  Koh,  Nk^NFFT,  NCP] =dvbtf rame2k; 
Kpilots=Koh ( : ,  l:Nk(2)); 

Xpilots=Xoh ( : ,  l:Nk(2)); 

[KH,  XpH,  IK,  Minv] =channelestimator (Kpilots ,  Xpilots , NFFT) ; 


%noise  variance  (account  for  BOOST  pilot) 

Ps  =  1;  %signal  power=unity  power=OFDM 

data  power 

L  =  NFFT/32;  %prefix  length (420) 

Ndata  =  1512; 

Ntps  =  17; 

Npilot  =  176; 


N  =  NFFT;  %OFDM  data  length (2112) 

r_pilot  =  16/9;  %pilot  power  boost  factor 

G  =  1  +  (L/  (Ndata+Ntps+Npilot’^r_pilot )  )  ; 

%eps_noise  =  Ps’^G/ ( log2  (M) ’^SDnoise^2 )  ;  %noise  power  (assume 
zero-mean) 


%compute  the  normalize  gain  factor  K (account  for  BOOST  pilot) 
N  =  NFFT; 

Ndata  =  1512; 

Npilot  =  176; 

Ntps  =  17; 

K  =  N/sqrt  (Ndata+Npilot’^r_pilot+Ntps )  ;  %gain  factor  to 
normalize  total  power 

%%  Ratios 

a 

o 

SNRdB  =  [-20,-15,-10,-5,0,5,10]; 

SIRdB  =-25^ (ones (size (SNRdB) )) ; 

Attn  =  1/10;  %  Attenuation  for  received  signal 

Rx_gain  =  le2;  %  Gain  for  shat 

SDint  =  1; 

SDsig_array  =  (10.^(SIRdB./20)).’^SDint; 

SDnoise_array  =  ( 10  .  (  (SNRdB+20 )  .  /  20 )  )  .  ’^SDsig_array  ; 

for  j  =  1  :  length ( SNRdB) 

SDsig  =  SDsig_array ( : , j ) ; 

SDnoise  =  SDnoise_array ( : , j ) ; 

eps_noise  =  Ps’^G/ ( log2  (M) ’^SDnoise^2  )  ;  %noise  power  (assume 
zero-mean) 

for  i  =  l:iter 

hlest_array  ( i ,  :  )  =hl  +  rho’^hl .  ( randn  (size  (hi )  )  +sqrt  ( - 

1 )  randn  (size  (hi)  )  )  / sqrt  (2); 

h2est_array  ( i ,  :  )  =h2  +  rho’^h2  .  ( randn  (size  (h2  )  )  +sqrt  ( - 

1)  ’^randn  (size  (h2)  )  )  /sqrt  (2)  ; 
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[numl_array ( i , : ) , num2_array ( i , : ) , den_array ( i , : ) ]  = 
kalmaninitial  (hlest_array ( i , : ) , h2est_array ( i , : ) , SDnoise ) ; 
end 

hlest  =  sum (hlest_array) /iter ; 
h2est  =  sum (h2est_array) /iter  ; 
numl  =  sum (numl_array) /iter  ; 
num2  =  sum (num2_array) /iter  ; 
den  =  sum (den_array) /iter  ; 

sim (model ) 

SNR(j,:)  =  lO’^loglO  ( var  ( s ) /var  (w)  ) +2  0  ; 

SIR(j,:)  =  lO’^loglO  (var  ( s )  )  /  ((var(yl)+var(y2))/2); 

Pb_weak(j,:)  =  err_weak(end,l); 

Pb_intl(j,:)  =  err_intl (end, 1 ) ; 

Pb_int2 ( j , : )  =  err_int2 (end, 1) ; 

toterr_weak ( : , j )  =  err_weak(size(err_weak,l)-lel:end,2); 

toterr_intl ( : , j )  =  err_intl(size(err_intl,l)-lel:end,2); 

toterr_int2 ( : , j )  =  err_int2 ( size (err_int2 , 1 ) -lei : end, 2 ) ; 


end 

RHO  =  rho  . ’^ones  ( size  ( SNR)  )  ; 

result  =  [RHO, SNR, SIR, Pb_intl , Pb_int2 , Pb_weak] ; 
str  =  'Simulation  complted  on  PC-2 ' ; 
email (str) ; 

2.  DVB-T2  FRAME  (courtesy  Prof.  Cristi) 

function  [Kmod,  Kdem,  Kdata,  Xoh,  Koh,  Nk,NFFT,  NCP] =dvbtf rame2k 
%  [Kmod,  Kdem,  Kdata,  Xoh,  Koh,  Nk,NFFT,  NCP] =dvbtf rame2k 
%  where 

o, 

o 

%  Kdata  =  data  FFT  indices  for  each  OFDM  frame  ("one"  based), 
68x1512 

%  Kmod  =  overhead  FFT  indices  for  each  OFDM  frame  ("one"  based), 
68x536 

%  with  Kmod= [Kdata,  Kpilots,  Ktps,  Knulls]  of  dimensions  68 

X  1512,  176,  17,  343 

%  Xoh  =  transmitted  values  of  each  overhead  subcarriers,  68x536 
%  with  Xoh= [Xpilots ,  Xtps,  Xnullls]  of  deimensions  68  x 

176,  17,  343 

%  Kdem  =  map  FFT  subcarriers  back  to  [data,  pilots,  tps,  nulls] 
vector 

%  Nk=[Ndata,  Npilots,  Ntps,  Nnulls] 


ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

oooooooooo 

%  From  the  DVB-T  Standard  (just  choose  one  mode  for  simplicity) 

%  2k  MODE  PARAMETRS: 

NFFT=2048; 
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NCP=NFFT/32;  %  length  of  Guard  Interval  (or  Cyclic  Prefix) 
Kmin=0;  Kmax=1704; 

Ndata=1512;  %  data  subcarriers 

Npilots=176;  %  pilots  subcarriers; 

Nnulls=343;  %  nulls  subcarriers 

Ntps=17;  %  Transmission  Parameter  Signalling  (TPS) 

subcarriers 

Nk=[Ndata,  Npilots,  Ntps,  Nnulls]; 


%  CELLS  ALLOCATION: 

%  Continual  Pilots  Cells 


kcp^ 

=  [0, 

LO 

00 

,  87,  141, 

156 

,  192 

r 

.  . 

201, 

255, 

279,  282, 

333, 

432, 

450,  ... 

483, 

525, 

531,  618, 

636, 

714, 

75 

i9,  ... 

765, 

780, 

804,  873, 

888, 

918, 

939,  ... 

942, 

969, 

984,  1050, 

1101,  11 

07, 

1110,  . 

1137 

,  1140 

,  1146,  12 

06, 

1269, 

1323,  1377 

1491 

,  1683 

,  1704]; 

%  Scattered  Pi 

lots  Cells 

Ksp 

(1:4: 

68,  :)  = 

ones (17,1) 

^  (3^ 

0  +  12^ 

(1: 

142) )  ; 

Ksp 

(2:4: 

68,  :)  = 

ones (17,1) 

^  (3^ 

1  +  12^ 

(0: 

141) )  ; 

Ksp 

(3:4: 

68,  :)  = 

ones (17,1) 

^  (3^ 

2  +  12’^ 

(0: 

141) )  ; 

Ksp 

(4:4: 

68,  :)  = 

ones (17,1) 

^  (3^ 

3  +  12^^ 

(0: 

141) )  ; 

%  ell  mod  4  =0 
%  ell  mod  4  =1 
%  ell  mod  4  =2 
%  ell  mod  4  =3 


%  TPS  (Transmission  Parameters  Signalling)  Cells 
ktps=[34,  50,  209,  346,  413,  ... 

569,  595,  688,  790,  901,  .  .  . 

1073,  1219,  1262,  1286,  1469,  ... 

1594,  1687]; 


%  Pseudoandom  BPSK  values  for  pilots,  depneding  on  subcarriers 
only) 

%  (for  simplicity  just  put  random  values) : 

Xp=  (4/3) ’^sign  (randn  (1,  NFFT)  )  ;  %  same  form  as  the  DVB-T  standard 

%  Xp=  (4/3)  ’^sqrt  (2)  ^  (sign  (randn  (1,NFFT)  )  +  sign  (randn  (1,NFFT)  ’^li)  )  ; 
%  Xp=  (0  :NFFT-1)  ’^11; 

%  Xp=ones (1,NFFT) ; 

%  Xp=l:NFFT; 


ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

g,  o,  g,  o,  g,  g, 
o  o  o  o  o  o 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

g,  g,  g,  g, 
o  o  o  o 

%  Translate  to  FFT  indices  (one  based) : 

%  Continual  Pilots 

kcp=kcp- (Kmax-Kmin) /2 ;  %  translate  to  FFT  indices  (0  based) 

I=f ind ( kcp<0 ) ; 
kcp (I) =NFFT+kcp (I) ; 

Kcp=ones  (68,1)  ’^kcp; 

%  Scattered  Pilots 

Ksp=Ksp- (Kmax-Kmin) /2 ;  %  translate  to  FFT  indices  (0  based) 

I=f ind (Ksp<0 ) ; 

Ksp  (I) =NFFT+Ksp (I)  ; 
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%  TPS  Carriers 

ktps=ktps- (Kmax-Kmin) /2 ;  %  translate  to  FFT  indices  (0  based) 

I=f ind ( ktps<0 ) ; 

ktps ( I ) =NFFT+ktps ( I ) ; 

Ktps=ones  (68, 1)  ’^ktps; 

%  Nulls 

kn=Kmax/2+l : (NFFT-Kmax/2  -1);  %  translate  to  FFT  indices  (0 

based) 

Kn=one  s  (68, 1 )  kn  ; 


%  Put  all  cells  together  in  a  matrix  Sk  68  x  NFFT : 

Sk=ones(68,  NFFT) ; 
for  ell=l:68 

Sk(ell,  Kn (ell , : ) +1 ) =0 ;  %  mark  nulls  with  0 

Sk(ell,  Kcp (ell , : ) +1 ) =-l ;  %  mark  pilots  with  -1 

Sk(ell,  Ksp (ell , : ) +1 ) =-l ;  %  same 

Sk(ell,  Ktps (ell ,:) +1 ) =-2 ;  %  TPS  carriers  with  -2 

end 

%  FFT  (one  based)  indices  for  data  (Sdata) ,  pilots  (Spilots) ,  TPS 
(Stps)  , 

%  nulls  (Snulls 
for  ell=l:68 

Idata=f ind (Sk (ell, :)==!) ;  Kdata (ell, : ) =Idata; 

Ipilots=f ind (Sk (ell, : ) ==-l) ;  Kpilots(ell,:)=Ipilots; 

Itps=find (Sk (ell,  : ) ==-2)  ;  Ktps (ell,  :)=Itps; 

Inulls=f ind ( Sk (ell , : ) ==0 ) ;  Knulls (ell, :)=Inulls; 

Xpilots (ell , : ) =Xp ( Ipilots ) ;  %  transmitted  pilots 

Xtps (ell , : ) =zeros ( 1 , Ntps ) ;  %  transmitted  TPS  (not  activated) 

end 

Xnulls=zeros ( 68 ,  Nnulls) ; 

%  Subcarriers  Overhhead 
Koh= [Kpilots ,  Ktps,  Knulls]; 

Xoh= [Xpilots ,  Xtps,  Xnulls] ; 

%  Modulator  mapping  from 

%  [data (68, 1512) ,  pilots ( 68 , 17 6 )  ,  tps(68,  17),  nulls(68,  343)] 

%  to  FFT  frame  subcarriers  (one  based) 

Kmod=[Kdata,  Kpilots,  Ktps,  Knulls]; 
for  i=l : 68 

[Ktemp,  Ik] =sort (Kmod ( i , : ) ) ; 

Kmod ( i , : ) =Ik; 

end 

%  Demodulator  mapping  from  FFT  subcarriers  back  to 
%  [data,  pilots,  tps,  nulls] 
for  ell=l : 68 

[xtemp,  I ] =sort (Kmod (ell , : ) ) ; 

Kdem (ell , : ) =I ; 

end 
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%  end 

3.  Channel  Estimation  (courtesy  Prof.  Cristi) 

function  [KH,  XpH,  IK,  Minv] =channelestimator (Kpilots , 
Xpilots,NFFT) 

%  [KH,  XpH,  IK,  Minv] =channeIestimator (Kpilots ,  Xpilots,  NFFT) 

%  To  test: 

%  clear 

%  [Kmod,  Kdem,  Kdata,  Xoh,  Koh,  Nk,NFFT,  NCP] =dvbtf rame2k; 

%  KpiIots=Koh ( : ,  l:Nk(2)); 

%  XpiIots=Xoh ( : ,  l:Nk(2)); 

%  I.  Parameters  for  Channel  Freq.  Response  Estimation  by  Pilot 
Response  Interpolation 
for  eii=I : 4 

K0=[Kpiiots (eii,  : ) ,  Kpiiots (eii  +  I,  :  )  , 

Kpiiots (eii  +  2,  : ) , Kpiiots (eii  +  3,  :  )  ]  ; 

X0=[ Xpilots (eii,  : ) ,  Xpilots (eii  +  I,  : )  , 

Xpilots (eii+2, : ) , Xpilots (eii+3, : ) ] ; 

[Ktemp,  I ] =sort (KO ) ; 

%  eliminate  repeated  subcarriers  due  to  Continual  Pilots 
Ki=K0(I);  Ki= [Ki,NFFT] ;  %  have  to  add  the  last 

subcarrier  NFFT 

dKi=dif f (Ki ) ;  Ii=f ind (dKi~=0 ) ;  %  find  non  repeated 

subcarriers 

KH(eii, : )=I (li) ; 

XpH (eii, : ) =X0 (I (li) ) ; 

end 

%  reorder  subcarriers  after  interpolation 
Ki=i:3:853; 

K2=2:3:85i; 

K3=3 : 3 : 852; 

K4=1197:3:2046; 

K5=ii98:3:2044; 

K6=ii99:3:2045; 

K7=2047; 

K8=2048; 

K9=854:il96; 

K= [ Kl , K2 , K3 , K4 , K5 , K6 , K7 , K8 , K9 ] ; 

[Ktemp,  IK] =sort (K) ; 

%  2 .  Matrix  for  Channel  Impulse  Response  Estimation 
%  Ail  pilots  values  and  indices 
Xaiip=zeros (I, NFFT) ; 
for  eii=i:4 

Xaiip (Kpiiots (eii, : ) ) =Xpiiots (eii, : ) ; 

end 

Kaiip=[0:3:852,  1196 : 3 : 2045] +1; 
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%  All  pilots  in  time  domain 
xallp=if ft (Xallp) ; 

%  Form  matrix 

L=64;  %  max  length  of  channel 

X=toeplitz (xallp,  [xallp(l),  xallp(NFFT:-l:2)]); 

X=X(:,1:L); 

F2048=if f t (eye (2048 )  )  ;  %  ifft  matrix 

F0=F2048 ( : ,  Kallp) ; 

[U,S,V]=svd(X) ; 

L0=55;  %  number  of  nonzero  singular  values  we  keep 

Minv=V(  :  ,  1  :L0)  ’^diag  (1 .  /diag  (S  (1  :L0,  1 :  LO )  )  )  (  :  ,  1  :L0)  ”^F0; 

%  end 

4.  Kalman  Filter  (courtesy  Prof.  Cristi) 

function  [numl , num2 , den]  =  kalmaninitial (hi , h2 , SDnoise ) 
coder . extrinsic ( ' ss2 tf '  )  ; 
coder . extrinsic ( ' dlqe '  )  ; 

%  hi  =  hi  /  hi (1) ; 
hi  =  [hi  0] ; 
h2  =  [0  h2] ; 

Nl=length (hi ) ; 

N2=length (h2 ) ; 

N=max ( [N1 , N2 ] ) -1 ;  %  system  order 

%  Form  state  space  matrices 
A=[-hl (2:N1) /hi  (1)  ; 

eye(N-l),  zeros (N-1 , 1 )] ; 

B=[l/hl(l) ; 

zeros (N-1 , 1 ) ] ; 

Bw=B  ; 

C=h2 (2 :N2) ; 

%  D  =  0; 

Q= (SDnoise/hl  ( 1 ) ) ^2 ; 

R=SDnoise"^2  ; 

%  Q  =  le-3; 

%  R  =  le-3; 

M  =  complex(zeros(length(A),l)); 

[M,P,Z,E]=dlqe (A,Bw,C,Q,R)  ; 

K=A’^M; 


Cx= [zeros ( 1 , N-1 ) , 1 ] ; 

N) 

%  pick  one 

of 

the 

states , 

say 

X  (n- 

%Cx=[l,  zeros ( 1 , N-1 )] ; 

%  pick  one 

of 

the 

states , 

say 

x(n) 

numl  =  complex ( zeros ( length (A) +1 , 1 ))  ; 
num2  =  complex ( zeros ( length (A) +1 , 1 )) ; 
deni  =  complex ( zeros ( length (A) +1 , 1 )) ; 

[numl,  deni  ]  =ss2  tf  (A-K’^C,  B,  Cx,  0,1); 
[num2,  den2  ]  =ss2 tf  (A-K’^C,  K,  Cx,  0,1); 
den  =  deni; 
end 
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